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ABSTRACT 

We have obtained /-band photometry of the neutron star X-ray transient Aql X-l 
during quiescence. We find a periodicity at 2.487 c d , which we interpret as twice 
the orbital frequency (19.30±0.05 h). Folding the data on the orbital period, we model 
the light curve variations as the ellipsoidal modulation of the secondary star. We 
determine the binary inclination to be 20°-31° (90 per cent confidence) and also 95 
per cent upper limits to the radial velocity semi-amplitude and rotational broadening 
of the secondary star to be 117kms _1 and 50kms _1 respectively. 

Key words: binaries: close - stars: individual: Aquila X-l - stars: neutron - X-rays: 
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1 INTRODUCTION 

Aquila X-1(=V1333 Aquilae) is a soft X-ray transient source 
that shows type I X-ray bursts (Koyama et al. 1981; Czerny, 
Czerny & Grindlay 1987), thereby indicating that the com- 
pact object is a neutron star. From quiescent observations 
the companion star has been identified to be a V=19.2 Kliv 
star (see Shahbaz, Casares & Charles and references within). 
Aql X-l is known to undergo regular X-ray and optical out- 
bursts on a timescale of ~1 year (Kaluzienski et al. 1977; 
Priedhorsky & Terrell 1984; Charles et al. 1980) much more 
frequently than the other neutron star transient Cen X-4 
(McClintock & Remillard 1990). 

Attempts to find the orbital period have revealed many 
modulations. Watson (1976) reported an unconfirmed 1.3 
day X-ray periodicity during the 1975 outburst. Chevalier 
& Ilovaisky (1991) have obtained an 18.97 hr periodicity 
from optical photometry during its active state, which they 
interpret as being the orbital period. 

Recently the RXTE All Sky Monitor showed Aql X-l 
to have undergone an X-ray outburst between late January 
and early March 1997 (Levine & Thomas 1997). Ilovaisky 
& Chevalier (1997) reported that Aql X-l was optically in 
quiescence by 30 March 1997. In this letter we report on 
our /-band photometry of Aql X-l obtained in June 1997, 
when the source was in quiescence. We are able to confirm 
the periodicity detected by Chevalier & Ilovaisky (1991) and 
determine the binary inclination. 

2 OBSERVATIONS AND DATA REDUCTION 

We monitored Aql X-l on the nights of 1997 June 20, 23, 
24, 25, and 26 UT, using the 1.3 m McGraw-Hill telescope 



at Michigan-Dartmouth-MIT (MDM) Observatory and a di- 
rect CCD camera equipped with a 1024 2 thinned Tektronix 
chip which yielded 0".51 pixeP 1 at f/7.5. The nights of 
June 21 and 22 were clear, but the source was too close to 
the bright moon to observe. To suppress scattered moonlight 
and to obtain maximum sensitivity to the cool secondary 
star, an interference filter approximating the Kron-Cousins 
/-band was used for all the exposures, which were generally 
480 s long. Except for a very few thin clouds on June 20, 
conditions were entirely photometric. The seeing was vari- 
able but generally usable (FWHM < 2"). Occasional focus 
changes and tracking problems spoiled a few frames. In or- 
der to obtain maximum leverage in the time series analysis 
we observed Aql X-l for as long as we could each night. 

For flatfield division we used median-filtered exposures 
of the twilight sky taken the same nights as the data. The 
flatfield pictures agreed very well from night to night, except 
for occasional ~ 1 per cent changes in the shadows cast by 
dust particles on the CCD window. Bias subtraction was 
accomplished using overscan regions of the picture. For all 
reductions we used IRAF^ routines. We measured eight stars 
(Table 1) in each picture using the apphot package; of a 
variety of software apertures, a l."4 radius gave the lowest 
scatter among the differential magnitudes, the averages of 
which are listed in Table 1. Photometric standardization was 
not attempted. The coordinates in Table 1 are derived from 
a fit of eight stars from the HST Guide Star Catalog vl.2 
in a single good-seeing frame, and are estimated accurate to 
< 0".3 from the scatter of the fit. The mean differential / 



* IRAF, the Image Reduction and Analysis Facility, is produced 
by the National Optical Astronomy Observatories. 
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Table 1. Differential magnitude and position of local field stars. 



RA 


JJcL 


AI 


C o rn in c lit s 


J ZUUU ) 


J ZUUU ) 


(mags) 




19 11 16.06 


+0 35 05.6 


2.272±0.022 


Aql X-l 


19 11 16.11 


+0 35 13.9 


(±0.108) 


local standard 


19 11 15.61 


+0 35 01.0 


2.002±0.018 


comparison star 


19 11 13.99 


+0 34 39.9 


-0.524±0.009 




19 11 17.03 


+0 35 30.5 


0.752±0.009 




19 11 17.00 


+0 35 22.0 


2.576±0.025 




19 11 18.63 


+0 35 24.0 


3.024±0.038 




19 11 16.83 


+0 34 55.9 


2.697±0.040 





magnitudes were computed by averaging over all frames and 
rejecting points more than three standard deviations from 
the median. For the local standard (~ 8" north of Aql X-l), 
only the standard deviation of the instrumental magnitudes 
is given. 



3 THE SEARCH FOR THE ORBITAL PERIOD 

We expect the quiescent optical modulation in Aql X-l to 
arise primarily from the ellipsoidal variations of the sec- 
ondary star, as it does in other SXTs. These variations arise 
since the observer sees differing aspects of the gravitationally 
distorted star as it orbits the compact object (van Paradijs 
& McClintock 1995). In theory the modulation should have 
two maxima and minima. The two minima may be unequal 
depending on the binary inclination, but the maxima should 
be equal. However, in practice, the light from the accretion 
disc (especially in the optical) contaminates the ellipsoidal 
variations of the secondary star, making detailed interpreta- 
tions of the optical light curves difficult (see Shahbaz, Naylor 
& Charles 1993). 

Therefore, we first analysed the optical light curve of 
Aql X-l using the phase dispersion minimization algorithm 
(Stellingwerf 1975). This technique is insensitive to the 
shape of the modulation, but does not remove the effects of 
the window function. The method groups the data in phase 
bins and seeks to minimize the dispersion within the bins. 
The deepest minimum of the statistic is the best estimate 
of the period. The PDM spectrum was computed in the fre- 
quency range 0.1 to 4.0 c d _1 at a resolution of 0.01c d _1 
with 20 phase bins. 

The common method of computing a discrete Fourier 
transform (DFT) and halving the estimate of the period 
is equivalent to assuming that the maxima and/or minima 
are of equal depths. This may not be the case, as the ob- 
served ellipsoidal modulation may contain unequal maxima 
and/or minima depending on the binary inclination and the 
contamination by the accretion disc. A Lomb-Scargle peri- 
odogram (Lomb 1976; Scargle 1982) of the data set was then 
computed with the same resolution and frequency range as 
was used for the PDM periodogram. 

Figure 1 shows the PDM and Lomb-Scargle period- 
grams, with three frequencies present at 0.524 c d~ x , 1.509 
c d _1 and 2.487 c d _1 . The suggested orbital period of 
Chevalier & Ilovaisky (1991) is marked at 18.97±0.02 hr. If 
this is the orbital period, then it should appear in the Lomb- 
Scargle periodogram as a peak at twice the frequency, i.e. 



at 2.530 c d . In both the PDM and Lomb-Scargle peri- 
odograms a periodicity at 2.487 c d _1 in present. Chevalier 
& Ilovaisky's period should also appear in the PDM peri- 
odogram, since the method is insensitive to the shape of the 
modulation in the light curve; a periodicity at 1.242 c d _1 
is present. We note that no periodicity at 1.242 c d _1 or 
2.487 c d _1 is present in the light curves of the comparison 
stars. We therefore assume the peak in the Lomb-Scargle 
periodogram at 2.487 c d _1 with a power of 31 to be real. 

In Figure 1 we also show the 99 per cent confidence 
level which allows us to demonstrate the significance on the 
peaks detected in the Lomb-Scargle periodogram. The level 
was calculated from a Monte Carlo simulation, which calcu- 
lates the maximum power in 10 000 sets of Gaussian noise 
with mean and variance equal to that of Aql X-l. Random 
peaks reaching a power of 7.8 are only found in 1 per cent 
of the artificial data sets, hence this defines our 99 per cent 
confidence level. Thus we conclude that the power level of 
the 2.487 c d _1 peak is substantially greater than 99 per cent 
level significant. In order to estimate the uncertainty in the 
2.487 c d" 1 peak position we created an artificial data set. 
with a mean, variance, semi-amplitude, amplitude difference 
between the minima and baseline equal to that in Aql X-l. 
The observed ellipsoidal modulation is well represented by 
a sinusoid of frequency 2/ which lags 90 degrees in phase, 
relative to sinusoid at frequency /. The amplitude of the si- 
nusoid at frequency / determines the difference between the 
minima in the final light curve. Since this difference is small 
(0.01 mags), the ellipsoidal light curve is very similar to a 
single sinusoid of frequency 2/. Using a Monte Carlo simu- 
lation we recorded the peak power in the Lomb-Scargle peri- 
odogram near the peak of interest, for the artificial data set. 
This was repeated 10 000 times to produce good statistics, 
giving a mean peak at 2.487 ±0.005c d~ (1-a uncertainty). 
Since this peak corresponds to twice the orbital frequency, 
the orbital period is then 1.244L0.003 c d _1 (=19.30±0.05 
h). 

Chevalier & Ilovaisky (1991) were taken over a 3 month 
baseline, when Aql X-l was in its active state. The V-band 
modulation they observe is probably due to a combination 
of X-ray heating of the secondary star, and of the accretion 
disc. Modulations of this kind are probably not coherent, 
due to the highly variable nature of the accretion disc, even 
more so in the V-band (van Paradijs & McClintock 1995). 
Our observations of Aql X-l were taken when the source 
was in quiescence, and where in the /-band the accretion 
disc contribution is negligible (Shahbaz, Casares & Charles 
1997). Therefore, what we see in the /-band is the charac- 
teristic double humped modulation per orbital cycle, due to 
the ellipsoidal shape of the secondary star. We therefore take 
the orbital period to be 19.30±0.05 h (=0.804±0.002 d). 

From Table 1 it can be seen that the scatter in the Aql 
X-l data is similar to stars of comparable brightness. This 
may suggest that the modulation we have detected is not 
real. However, it should be noted that we can still use the 
largest amplitude of an ellipsoidal modulation that could 
conceivably be hidden within the noise in the data, to con- 
strain the inclination of the system. 
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4 MODEL FITTING AND DERIVED 
PARAMETERS 

Using an orbital period P=19.30 h (=0.804 d) we folded the 
Aql X-l and comparison star data. As there is no spectro- 
scopic ephemeris defining phase 0.0 (i.e. superior conjunc- 
tion of the secondary star) we used an arbitrary value for 
To. We then shifted the light in phase so that the deeper 
minima corresponded to phase 0.0. [Note that the choice 
is arbitrary. In the fitting procedure (see later), one of the 
free parameters is a phase shift.] The data were then binned 
into 15 phase bins, and the standard error and mean were 
then calculated for each bin after discarding the data points 
which were more than 3-cr away from the mean of each bin. 
This process removed 4 and 12 discrepant data points from 
the Aql X-l and comparison star data respectively; the total 
number of data points was 152. Figure 2 shows the resulting 
light curves for Aql X-l and a comparison star. 

We fitted the orbital /-band light curve of Aql X-l with 
an ellipsoidal model, similar to that used for the other ob- 
jects in our programme (see Shahbaz, Naylor & Charles 1997 
and references within). The model describes the light curve 
generated by a Roche lobe filling star, where each element of 
area on the surface of the star is assumed to emit blackbody 
radiation. The temperature distribution over its surface is 
assumed to vary according to Von Zeipel's (1924) gravity 
darkening law and a linear limb darkening law is used. For 
a detailed description of the model see Shahbaz, Naylor & 
Charles (1993). The model parameters were the binary mass 
ratio (q —M1/M2 where Mi and M2 are the mass of the 
compact object and secondary star respectively), the incli- 
nation (i) and the effective temperature of the secondary 
star (T e ft). A bright spot, i.e. the region where the accretion 
stream hits the accretion disc was also included. The free 
parameters of the model were a phase shift, the normalisa- 
tion of the bright spot and the normalisation of the light 
curve. Thus the fitting procedure is not dependent on the 
choice for To. We performed least-squares fits to the data 
using this model, grid searching q in the range 1-10 and i in 
the range 5°-50°. A T cff of 4620 K, appropriate for a KllV 
star was used. We used a gravity darkening exponent of 0.05; 
Sarna (1989) finds the gravity darkening exponent for con- 
vective stars with M2 <0.7 to be 0.05. The limb darkening 
coefficient for the J-band (9000A) and the varying temper- 
ature around the surface of the secondary star was interpo- 
lated from the tables given by Al-Naimiy (1978). The best 
fit gives a Xv °f 1-45 at q = 3, i = 25°, with a phase shift 
of 0.0033</>(=0.027 d). This fit is plotted in Figure 2 as the 
dashed line. The phase shift allows us to estimate To to be 
HJD 2450623.450±0.016. 

Figure 3 shows the resulting (i, q) diagram with the 
68 and 90 per cent confidence regions marked, calculated 
according to Lampton, Margon & Bowyer (1976) for three 
parameters (the phase shift, bright spot normalisation and 
the light curve normalisation), after the error bars had been 
rescaled to give a fit with a of 1- As one can see, the 
ellipsoidal variation is only weakly dependent on q. The 68 
and 90 per cent confidence regions for the binary inclination 
are 21°-30° and 20°-31° respectively. 

Shahbaz, Casares & Charles (1997) determine the frac- 
tion of light arising from the accretion disc to be 6±3 per 
cent at 6000A. With increasing wavelength the disc contam- 



ination decreases (see McClintock & Remillard 1986; Shah- 
baz, et al. 1996 and references within), therefore the /-band 
disc contamination should be lower. If we assume a 95 per 
cent upper limit of 9 per cent for the disc contamination, 
the effect of this on the binary inclination is to increase it 
by at most 1°. 

From Figure 2 it can be seen that there seems to be 
some systematic trend in the data of the comparison star. 
As noted earlier the scatter in the comparison star data is 
comparable to that in the Aql X-l data (see also Table 1). 
If we assume that the modulation we observe in the Aql 
X-l data is not due to the ellipsoidal modulation of the 
secondary star, and what we have fitted is the maximum 
possible modulation that could be present, then our fits al- 
low us to determine a 95 percent upper limit of 31° to the 
binary inclination. 



5 DISCUSSION 

5.1 The period difference 

The quiescent period (19.30 h) we find is slightly different 
(~ 2 per cent) from the outburst period (18.97 h) obtained 
by Chevalier & Ilovaisky. Outburst and quiescent period dif- 
ferences are observed in SXTs (see O'Donoghue & Charles 
1996), but the outburst modulation (caused by the preces- 
sion of the accretion disc; Whitehurst 1988) is usually a 
few per cent longer than the orbital period. If we inter- 
pret Chevalier & Ilovaisky's outburst period as being due 
to the superhump phenomenon, then the superhump pe- 
riod is shorter than the orbital period by 2 per cent, i.e. 
we have a "negative" superhump. This may seem unusual, 
but negative superhumps have previously been observed in 
cataclysmic variables, e.g. in V503 Cyg the superhump pe- 
riod is shorter than the orbital period by 3 per cent (Harvey, 
Skillman, Patterson & Ringwald 1995). 



5.2 Upper limits on K2 and «sin i 

We can compare Aql X-l to the other long period neutron 
star SXT Cen X-4. Assuming that our inclination is correct, 
and that both systems have similar mass neutron stars, then 
one can estimate K2, the semi-amplitude of the radial veloc- 
ity curve for Aql X-l (K2 scales with sin i). Using 7^2=146 
km s _1 (McClintock & Remillard 1990) and i = 40° for Cen 
X-4 and i < 31° for Aql X-l, we find K 2 < 117kms _1 (95 
per cent upper limit) for Aql X-l 

We can also determine the upper limit to the rota- 
tional broadening of the secondary star (v sin i) assuming 
that the secondary star fills its Roche lobe and its spin is 
tidally locked to the binary period. Eliminating K2 from the 
mass function and the equation which relates the rotational 
broadening to the radial velocity semi-amplitude and mass 
ratio (Wade & Home 1986), we obtain the expression 

t>sin i = 283sin i{M 2 /P) 1/3 km s _1 . (1) 

The secondary star must be less massive than a main se- 
quence star of the same spectral type, 0.8 M@. Then using 
equation (1) with P=19.3 hr and i < 31° gives «sin i < 
50 km s -1 (95 per cent upper limit). This value for i>sin i is 
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consistent with the limits obtained by Shahbaz, Casares & 
Charles (1997). 



6 CONCLUSION 

Using /-band photometry of Aql X— 1 taken during quies- 
cence, we have detected a periodicity at 2.487 c d~ x , which 
we interpret as twice the orbital frequency. Folding the data 
on the orbital period (19.30±0.05 hr), we fit the ellipsoidal 
modulation of the secondary star. We determine 95 per cent 
upper limits to the binary inclination, radial velocity semi- 
amplitude and rotational broadening of the secondary star 
to be 31°, 117kms _1 and 50 km s _1 respectively. 
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Figure 1. Results of the period search using the /-band data of Aql X-l during quiescence. Top: Phase Dispersion Minimization (PDM) 
periodogram. The frequency 1.242 c d~ 1 which we interpret as the orbital frequency is marked with a star; twice this frequency is also 
shown (white circles). Middle : Window function. Bottom : Lomb-Scargle power spectrum. The vertical dashed line is the orbital period 
determined by Chevalier & Ilovaisky (1991) at 18.97 hr. A Monte Carlo simulation provides the 99 per cent confidence level, shown as 
the dashed horizontal line. We interpret the periodicity at 2.487 c d -1 , present in the Lomb-Scarglc and PDM pcriodograms, as twice 
the orbital frequency (see text), since it is consistent with twice Chevalier & Ilovaisky 's orbital frequency. 
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Figure 2. The phase-binned light curves for Aql X— 1 (top) and a comparison star (bottom). In the top panel we also show the phase- 
folded data for Aql X-l (light weighted crosses). The dotted line in the comparison star light curve is a linear fit. The best fit ellipsoidal 
model (dashed line) at q = 3 and i = 25° is shown in the light curve of Aql X-l. We show 1.5 orbital cycles for clarity. 
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Figure 3. The 68 and 90 per cent confidence level solutions for the ellipsoidal model fits to the Aql X-l /-band light curve. The 
marks the best solution. 
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